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High-fidelity gate operations for quantum computing beyond dephasing time limits. 
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The implementation of quantum gates with fidelities that exceed the threshold for reliable quan¬ 
tum computing requires robust gates whose performance is not limited by the precision of the 
available control fields. The performance of these gates also should not be affected by the noisy 
environment of the quantum register. Here we use randomized benchmarking of quantum gate oper¬ 
ations to compare the performance of different families of gates that compensate errors in the control 
field amplitudes and decouple the system from the environmental noise. We obtain average fidelities 
of up to 99.8%, which exceeds the threshold value for some quantum error correction schemes as 
well as the expected limit from the dephasing induced by the environment. 
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Scalable quantum computing requires gate operations 
with fidelities above a certain threshold [l|, . Reaching 

this threshold remains challenging, primarily due to two 
issues: (i) the control fields driving the gate operations 
have often experimental uncertainties, such as amplitude 
errors, that are bigger than the allowed deviations and 
(ii) the quantum systems cannot be completely shielded 
from the effects of a noisy environment. Different tech¬ 
niques have been developed to overcome these obstacles, 
such as dynamical decoupling (DD) for suppressing the 
effect of environmental noise and thus extend the coher¬ 
ence time of the system. An ideal dynamical decoupling 
sequence completely eliminates the interactions of the 
system with its environment and thereby “freezes” the 
system, apart from the evolution under the internal sys¬ 
tem Hamiltonian. This approach can therefore imple¬ 
ment high-fidelity quantum memories. If, however, the 
specific application requires that the system evolves un¬ 
der a suitable control Hamiltonian, such as in a quan¬ 
tum information processor, the protection scheme and 
the control operations must be applied simultaneously. It 
is then necessary to design gates that combine processing 
with decoupling in such a way that the effect of the con¬ 
trol fields survives, while all unwanted interactions are 
eliminated Q. 

Measuring the fidelity of gate operations in the region of 
the threshold for scalable quantum computation cannot 
be done on individual gates, since the measurement errors 
can be comparable to or greater than the gate errors. For 
this purpose, the “randomized benchmarking” procedure 
was developed It measures an average gate fidelity 
of a random selection of gates. This allows a much more 
precise determination of the average gate fidelity. It does 
not provide information about the fidelity of individual 
gates, but it provides estimates of the fidelity of sequences 
of gates, which is the relevant quantity for the analysis 
of large-scale computation. 

In this paper, we present several families of gate oper¬ 
ations that combine robustness against unwanted varia¬ 
tions in the amplitudes of the control fields (flip angle 


errors) with protection against environmental noise. For 
this purpose, we use Clifford gates, which allow universal 
quantum computation when combined with magic states 
[5l-l7|. We determine their fidelities experimentally, us¬ 
ing a nuclear spins as qubit and a different nuclear spin 
system as the noisy environment. 

We consider a system qubit S, which is coupled to a noisy 
environment via a dephasing interaction 

'Use = b{t)Sz, (1) 

where the (semi-)classical field b{t) has a finite correlation 
time Tc- Dynamical decoupling by sequences of inversion 
pulses can suppress the dephasing effect of this interac¬ 
tion and prevent the decay of the coherence in the system 
[1, Q . This is well explored in the context of quantum 
memories, where the goal is to keep a specific quantum 
state unchanged. Here the effect of DD can be roughly 
described as a prolongation of the dephasing time T 2 . 

In the case of quantum computing, the goal is not the 
preservation of a quantum state, but the precise control 
of a quantum system such that it evolves along a well- 
defined path in Hilbert space. Dephasing then creates 
errors roughly oc (1 — Ri t/T 2 , where r is the du¬ 

ration of the gate operation and T 2 the dephasing time. 
Extending the dephasing time by DD therefore helps also 
for quantum computing. However, in this case, an ad¬ 
ditional complication arises: the dynamical decoupling 
sequences that are used for protecting quantum memo¬ 
ries would equally decouple the system from the external 
fields that are applied to control the system’s path. It is 
therefore essential to combine dynamical decoupling and 
gate operations in such a way that they do not interfere 
destructively. Two general approaches for this have been 
proposed for solving this problem: t o ap ply gate opera¬ 
tions and DD operations successively [1^ , or to interleave 
gate operation and DD by splitting the gate operation 
into as many elements as there are delays between the 
DD pulses and modify them in such a way that the effect 
of the DD pulses is to revert these modifications and the 
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overall effect becomes that of the targeted gate operation 

A third option is to design gate operations that are ro¬ 
bust against experimental imperfections and have ‘built- 
in’ the effect of a dynamical decoupling sequence. Such 
a decoupling effect is, e.g., built into 7r-pulses around an 
axis in the a;?/-plane: they invert the system operator Sz 
and therefore the system environment interaction given 
in eq. ©• The sequences that we consider for random¬ 
ized benchmarking consist to 50% of 7r-pulses. Accord¬ 
ingly, they automatically reduce the effect of the system- 
environment interaction, although their efficiency may be 
lower than that of specialized DD-sequences. 

Apart from reducing the coherence time, high fidelity 
operations also need to be robust against errors in the 
control fields driving the gate operations. The BBl com¬ 
posite pulse (see figure [T^) was introduced by Wimperis 
[l^ as a possible scheme for generating composite rota¬ 
tion pulses that are well compensated against amplitude 
errors. The sequence for generating a compensated rota¬ 
tion by an angle 9 around an axis in the ly-plane is 

Rip {^9^ Ri^j^ip (^TT^ R^f^-^-ip (27r) (^) • 

Here, ip defines the orientation of the rotation axis and 
/3 = cos“^ (—0/47r). For the present purpose, we consider 
the cases 9 = tt/2 and tt, where (3go ~ 1.7 and /3i8o ~ 1-8. 
If the four 7r-pulses of this sequence are separated in time 
by a delay that is twice the duration of the initial d-pulse, 
this yields a sequence that is not only robust to amplitude 
errors, but also to environmental noise: it conforms to the 
“compute, then decouple” approach and the four 7r-pulse 
generate the DD cycle. 

A robust TT-pulse can also be generated by concatenating 
5 TT-pulses with the phase 7r/6, 0,7r/2, 0, tt/G [l^ (see fig¬ 
ure m). This generates an inversion of the z-component, 
but in addition also a —tt/S rotation around the z-axis. 
We will refer to this pulse as KDD-5 = Rz{—tt/3)Ro = 
RqRz{tt/3), where Rq is the targeted 7r-rotation around 
the x-axis. If it is used in a sequence like randomized 
benchmarking, the additional rotation can be taken 
into account by adjusting the coordinate system. 

To determine the average fidelity of the gate operations, 
we use the scheme proposed by [l^ [2^ . It requires the 
application of tt/2 and tt rotations around the axes of 
the coordinate system. The individual gates are 1-qubit 
Clifford gates C = PG, where P indicates a unit oper¬ 
ation or a TT-rotation around one of the coordinate axes 
(8 different operations) and G a. tt/2 rotation around a 
coordinate axis (6 different operations, 48 total). Rota¬ 
tions around ±z axes are implemented through a change 
to the rotating frame definition [l^, [2^ . At the end of 
the sequence a recovery operation is applied that corre¬ 
sponds to the inverse of the sequence up to this point, so 
that the full sequence becomes a unit operation in the ab¬ 
sence of imperfections. For optimal randomization, the 
recovery operation itself consists of two random P oper¬ 
ations sandwiching the complement i?, which is another 


Clifford gate. For a single qubit, the effect of a non-ideal 
identity operation can be written as 

. , d 

Pout ~ (1 d)pin “1. 

Here, pin is the initial state, 1 the unit operator and d is 
the depolarizing parameter. Here, we use the fact that for 
a sequence of Clifford gates, the result can be represented 
as a depolarization channel [U [22| . The trace fidelity is 
then 

d 

F = tr{poutPin\ = 1 “ 2' 
and the error per gate is the difference 

EPG =\-F =-. 
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If the error per gate is purely random, we expect that 
the average fidelity for sequences with m gate operations 
decays as 

F=i(l + (l-d)™)«i(l + e—''). 

For the experiment, we used the ^^C nuclear spins of 
adamantane, a molecular crystal. The protons in the 
same crystal provide the noisy environment, which leads 
to a dephasing time of 360 ps. A single refocusing pulse 
(Hahn echo) can extend the dephasing time to 740 ps. 
The ultimate limit on the information storage in this 
system is given by the energy relaxation time ri=1.52 
s. 

The system was initially in thermal equilibrium, with 
Pin oc Sz- For each measurement, we averaged over a 
set of 32 random sequences composed of I to 80 gate 
operations. After the restore operation, the signal was 
measured by applying one additional readout pulse that 
converted the remaining Sz component of the density op¬ 
erator into transverse magnetization. The free precession 
signal was acquired, Fourier-transformed, and the sig¬ 
nal corresponding to the CH 2 carbon was integrated and 
used as a measure for the survival probability 
Figure [1] gives an overview over the different sequences 
tested in this work. First (not shown in the figure) we 
tested gates without DD protection, using two different 
implementation for P and G rotations: i) simple rectan¬ 
gular pulses and ii) BBI pulses, which compensate errors 
in the amplitudes of the control fields. In the second ap¬ 
proach (H^) we generated a protected G rotation by sep¬ 
arating the four 7r-pulses of the BBI composite rotation 
by a delay that is twice the duration of the initial tt/2- 
pulse. The P rotation in this scheme was implemented 
by one KDD-5 pulse. The third set of sequences (Hb-e) 
was designed to reduce experimental imperfections and 
decouple the system from its environment by interleaing 
the P and G rotations with dynamical decoupling pulses 
in different ways, using BBI pulses for all PG rotations. 
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FIG. 1; Different sequences of gate operations compared in 
this work. The donble arrows indicate the gate dnration r, 
defined as the time period between the beginning of one P 
gate to the beginning of the next P gate, a) A protected G 
rotation is created by separating in time the four 7r-pulses 
of the BBl composite rotation by a delay that is twice the 
duration of the initial 7r/2-pulse, while the P rotation was 
implemented by one KDD-5 pulse. From b) to e) PG gates 
were implemented by BBl pulses and interleaved in different 
ways with the DD sequence XY-16 (see text), f) and g) show 
the KDD-5 and BBl pulses. 


In b) and c) the DD pulses were also implemented as 
BBl pulses, while in d) and e) the DD pulses were im¬ 
plemented by rectangular pulses. 



FIG. 2: Experimental survival probabilities for sequences of 
BB 1-gates as a function of the number of steps. The circles 
represent the survival probability averaged over 32 gate se¬ 
quences. The solid lines are represent an exponential fit and 
theoretical predictions (see text) . 

Figure [5] shows the experimental results obtained when 
composite BBl pulses were used to implement the Clif¬ 
ford gates. Here, the gate time was 76 /rs. The blue 
and green lines show the theoretical prediction for pure 
dephasing with times of 750 /iS and 340 /rs, respectively, 
which would result in EPGs of 3.2%, and 6.7%. The 
straight lines through the experimental data points are 
fits to exponential decays, the corresponding EPG is 
0.32 ± 0.03%, which is significantly better then the EPG 


FIG. 3: Experimental results for different type of gates: rect¬ 
angular pulses (red), BBl-gates (blue), a compute then de¬ 
couple approach built from BBl pulses and KDD-5 pulses 
(green), and protected gates by dynamical decoupling (black). 


extrapolated from the T 2 values. This is a clearly indica¬ 
tion that the environment is refocused during the bench¬ 
marking experiment, even without applying DD pulses. 
In Figure [3] we compare the BBl pulses with the three 
other cases: i) The approach defined in figure[T^, ii) inter¬ 
leaving BBl pulses with the DD sequence XY-16 (figure 
[TJ:) and iii) rectangular pulses. The gate time of the of 
the hrst scheme is 88 ^s and its performance is compara¬ 
ble to that of BBl, the observed EPG is 0.34 ±0.03%. In 
the case of interleaving BBl pulses with DD, the gate 
time is 152 /rs, which is approximately twice the gate 
time for BBl gates without DD protection. However the 
observed EPG is 0.22 ± 0.03%, showing that dynamical 
decoupling sequences can further increase the fidelity of 
quantum gates. The rectangular pulses do not exhibit 
an exponential decay, since those pulses are not built to 
compensate r.f. field inhomogeneities. In this case we 
cannot assign a single value to the EPG. For a quanti¬ 
tative evaluation, we fitted the experimental data to the 

k 

function o'" . The parameters obtained from this fit are 
a ~ 0.88 and k 0.44 

Figure H] compares the decay of the survival probability 
vs. time for the case of a benchmarking sequence whose 
gates are protected against environmental noise by an 
XY-4 DD sequence with the decay of coherence if only 
DD is applied, with no interleaved gate operations. For 
the benchmarking experiment the initial state was ther¬ 
mal equilibrium, with pin oc Sz, while in the pure DD 
experiment the system was initialized in a superposition 
state. Clearly, the benchmarking results in much slower 
decays then the XY-4 sequence. This result can be traced 
to the fact that the pulse imperfections in an XY-4 cy¬ 
cle result in a propagator Ui ~ 1 + Q, where Q is the 
error term. In the case of DD, the overall propagator is 
U4n = 114(1 + Q)” ~ 1 ± nQ : for small number of cycles 
n, the error grows linearly. In the case of benchmark¬ 
ing, the gates that are interleaved with the DD sequence 
modulate the error in a random manner, resulting in a 
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FIG. 4: Decay of the survival probability measured by inter¬ 
leaving quantum gates with the dynamical decoupling XY-4 
and pure dynamical decoupling, without gate operations. In 
the case of the XY-4 sequence, the decay can be fitted by 
a quadratic function (—0.02t^ — 0.64t -|- 0.97). The bench¬ 
marking and XY-16 results were fitted with the exponential 
functions and respectively 

random walk in Liouville space. As a result, an increase 
oc n may be expected for a random sequence. This is 
an expression of the general principle that in most cases, 
a decoupling scheme exists that is better than a simple 
repetition of the same basic cycle. 

The experimental data clearly show that a robust de¬ 
sign of the gate operations can make their performance 
largely independent of experimental imperfections. This 
does not mean that the error per gate can be made ar¬ 
bitrarily small, since not all interactions that contribute 
to the errors are actually suppressed. In particular, two 
types of environmental noise cannot be refocused by ex¬ 
plicit or implicit DD: (i) interactions that are not linear 
in the qubit operators S and (ii) interactions whose corre¬ 
lation time is shorter than the experimentally accessible 
time scales. In the present system, couplings 

are an example of (i): they are bilinear in the spin in¬ 
teractions and are not eliminated by the DD sequences 
used here, although other sequences exist for suppress¬ 
ing them [ 2 ^. Processes of category (ii) are those 
that contribute to the energy relaxation (Ti-processes): 
their correlation times must be comparable to or shorter 
than the Larmor precession period. Their contributions 
can be estimated from available experimental data. The 
longest dephasing times that have been reached in this 
system for isotropic decoupling sequences indicate that 
Ti 50 ms [ 2 ^. This value was achieved by using ro¬ 
bust pulses and interpulse delays close to zero and it is 
also in qualitative agreement with numerical estimates 
for the broadening by dipolar coupling in natural 

abundance adamantane. The contribution from fast pro¬ 
cesses can be estimated from the energy relaxation time 
as T 2 =Ti K, 1.52 s. We can therefore neglect processes 
of type (ii). This indicates that further improvements 
would require pulses that compensate also for homonu- 
clear (bilinear) interactions. If the remaining dephasing 


Gate 

BBl 

(a) 

(b) 

(c) 

(d) 

(e) 

r[ps] 

76 

88 

116 

152 

336 

384 

EPG^[10-4] 

5 

6 

8 

10 

22 

25 

EPGm[10-4] 

317 

364 

472 

604 

1191 

1322 

EPGexp [10-4] 

32±3 

34±3 

28±3 

22±3 

172±6 

47±3 


TABLE I: Summary of the gate operations tested. For each 
type of gate, the duration r for a PG operation is given, the 
limiting error per gates, EPGm and EPGm, which is expected 
for pure dephasing with a dephasing time measured by Hahn 
Echo and pure dynamical decoupling, respectively. The labels 
a-e refers to the gates types illustrated in the figure [T] . 


rate is r 2 this should contribute an error per gate of 
^r/Ta 


Table U summarizes the results. For each type of gate 
it shows the observed EPGexp- All gate operations in¬ 
clude decoupling elements and thus result in an error per 
gate that is lower than the error rate EPGm estimated 
from the dephasing time measured by Hahn Echo. This 
is a clear indication that gate operations alone provide 
some dynamical decoupling. Each sequence is also com¬ 
pared to the minimum EPGm values we could reach in 
the present system, EPGm ~ T/T 2 , which is expected 
from pure dephasing due to the remaining environmen¬ 
tal noise. These values were obtained by DD with ro¬ 
bust pulses and very short delays between the pulses 
. The best performance was achieved by the gate type 
(c), which corresponds to BBl gates protected by the DD 
sequence XY-16. 

The implementation of robust high-fidelity gate opera¬ 
tions is an essential step towards reliable and scalable 
quantum computing. In this work we use randomized 
benchmarking of single qubit quantum gates to compare 
the performance of different families of gates that com¬ 
pensate errors in the control field amplitudes and decou¬ 
ple the system from the environmental noise. In some 
cases the total duration of the experiments (from~ 1 ms 
to ^ 30 ms) exceeds the dephasing time, measured by 
Hahn echoes (760 ps), by almost two orders of magni¬ 
tude. We could obtain average fidelities exceeding the 
expected limit from the dephasing induced by the en¬ 
vironment. This is a clear indication that the effect of 
the noisy environment is reduced by the sequence of ap¬ 
plied gate operations. The best average fidelity observed, 
99.8% , was achieved by interleaving gate operations im¬ 
plemented by composite pulses with the dynamical de¬ 
coupling sequence XY-16. 

The improvements in gate accuracy by decoupling meth¬ 
ods, as observed in this work, implies a reduction of the 
overhead cost of QEG since more noise can be tolerated 
by QEG codes combined with DD than by QEG alone [l^. 
Therefore, future work will be devoted to investigate the 
performance of QEG codes combined with dynamical de- 
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coupling sequences and to test the performance of DD 
protected gates in systems with multiple qubits. 
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